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Abstract. Using the coupled cluster method we investigate spin-s J1-J2 Heisenberg antiferromagnets 
(HAFs) on an infinite, anisotropic, two-dimensional triangular lattice for the two cases where the spin 
quantum number s = 1 and s = |. With respect to an underlying square-lattice geometry the model has 
antiferromagnetic (Ji > 0) bonds between nearest neighbours and competing (J2 > 0) bonds between 
next-nearest neighbours across only one of the diagonals of each square plaquette, the same diagonal in 
each square. In a topologically equivalent triangular-lattice geometry, the model has two types of nearest- 
neighbour bonds: namely the J2 = kJi bonds along parallel chains and the Ji bonds producing an inter- 
chain coupling. The model thus interpolates between an isotropic HIAF on the square lattice at one limit 
(k = 0) and a set of decoupled chains at the other limit {k — >■ cxd), with the isotropic HAF on the triangular 
lattice in between at «: = 1. For both the spin-1 model and the spin-| model we find a second-order type 
of quantum phase transition at Kc = 0.615 ± 0.010 and Kc = 0.575 ± 0.005 respectively, between a Neel 
antiferromagnetic state and a helically ordered state. In both cases the ground-state energy E and its 
first derivative dE/dn are continuous at «: = Hc, while the order parameter for the transition (viz., the 
average ground-state on-site magnetization) does not go to zero there on either side of the transition. The 
phase transition at n — Kc between the Neel antiferromagnetic phase and the helical phase for both the 
s = 1 and s = I cases is analogous to that also observed in our previous work for the s = § case at 
a value Kc = 0.80 ± 0.01. However, for the higher spin values the transition appears to be of continuous 
(second-order) type, exactly as in the classical case, whereas for the s = | case it appears to be weakly 
first-order in nature (although a second-order transition could not be ruled out entirely). 

PACS. 75.10.Jm Quantized spin models - 75.30.Kz Magnetic phase boundaries - 75.50.Ee Antiferromag- 
netics 



1 Introduction 

, In recent years, the theoretical study of two-dimensional 
(2D) quantum spin systems has been intensely motivated 
\ by the fact that such models often describe well the prop- 

■ erties of real magnetic materials of great experimental in- 
' terest. It is an encouraging fact that experiments have of- 
, ten supported theoretical predictions or vice versa. More- 

■ over, the study of frustration and quantum fluctuation in 
\ quantum spin-lattice systems has developed into an ex- 
tremely active area of research. The interplay between 
frustration and quantum fluctuations in magnetic systems 
can produce a wide range of fascinating quantum phases 
[1112113] • Both effects are in principle capable of destabilis- 
ing or completely destroying the magnetic order of the spin 
system. In turn this might then lead to the formation of a 
spin-liquid phase or be responsible for other quantum phe- 
nomena of similar significant interest. The driving force 
for the differentiated manifestation of the various kinds of 
possible quantum effects can, in principle, also come from 
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the type and nature of the underlying crystallographic lat- 
tice, from the number and variety of the magnetic bonds, 
and from the magnitude of the spin quantum numbers of 
the atoms that reside on the atomic lattice sites |3]. It 
is thus of considerable interest to try to study the effects 
of each of these driving forces in turn for specific model 
Hamiltonians. 

A particularly well studied 2D model is the frustrated 
spin-i J1-J2 model on the square lattice with nearest- 
neighbour (NN) bonds (Ji) and next-nearest-neighbour 
(NNN) bonds (J2), for which it is now well accepted that 
in the case where both sorts of bonds are antiferromag- 
netic (Ji > 0;i = 1,2) there exist two antiferromagnetic 
phases exhibiting magnetic long-range order (LRO) at 
small and at large values of a = J2/J1 respectively. These 
are separated by an intermediate quantum paramagnetic 
phase without magnetic LRO in the parameter regime 
< a < , where ~ 0.4 and ctcj w 0.6. For 
a < the ground-state (gs) phase exhibits Neel mag- 
netic LRO, whereas for a > it exhibits collinear stripe- 
ordered LRO. 
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As already noted above, the spin quantum number s 
can, both in principle and in practice, play an important 
role in the phase structure of strongly correlated spin- 
lattice systems, which often display rich and interesting 
phase scenarios due to the interplay between the quantum 
fluctuations and the competing interactions. Varying the 
spin quantum number s can tune the strength of the quan- 
tum fluctuations and lead to fascinating phenomena 5^. 
A well-known example of such spin-dependent behaviour 
is the gapped Haldane phase [6] in s = 1 one-dimensional 
(ID) chains, which is not present in their s = ^ counter- 
parts. 

Some recent studies on large-spin (i.e., s > ^) systems 
include: (a) the comparison of the Heisenberg antiferro- 
magnet (HAF) on the Sierpihski gasket with the corre- 
sponding HAFs on various regular 2D lattices, including 
the square, honeycomb, triangular and kagome lattices, 
for the cases s = i, 1, and | [7]; (b) the one-dimensional 
(ID) Heisenberg chain for both integral and half-odd in- 
tegral values of s up to a value of 10 [51IU1ITU]: (c) the 2D 
J- J' model on the square lattice containing two different 
types of NN bonds, for values of s between ^ and 2 (Fj; 
d) the 2D spin-anisotropic XX Z Heisenberg Hamiltonian 
for s — 1 p fm[T2l[T3] : e) the spatially anisotropic Ji-Ji'- 
J2 model on the 2D square lattice for values of s up to 
4 [TillTS]: (f) the spin-anisotropic J^^^-J^^^ model for 
s = 1 [iB]; (g) the pure J1-J2 model for s = 1 [17]; (h) the 
2D Union Jack model for s = 1 and s = | JH]; and the 
2D Heisenberg model on the honeycomb lattice for s = 1 
[19j . Another example that has been experimentally stud- 
ied involves the investigation of the single-ion anisotropy 
energy in the 2D kagome lattice for the case s = | pp] . 

Also noteworthy in this context is the recent discov- 
ery of superconductivity with a transition temperature at 
Tc « 26 K in the layered iron-based compound LaOFeAs, 
when doped by partial substitution of the oxygen atoms by 
fluorine atoms [2T| . This finding has been followed by the 
rapid discovery of superconductivity at even higher values 
of Tc (> 50 K) in a broad class of similar quaternary com- 
pounds. Enormous interest has thereby been engendered 
in this class of materials. The very recent first-principles 
calculations [351 shows, for example, that the undoped 
parent precursor material LaOFeAs of the first material 
investigated in this oxypnictide class is well described by 
the spin-1 J1-J2 model on the square lattice. 

We have previously used the coupled cluster method 
(CCM) [23..24..25, to study the magnetic order in a spin- 
half interpolating square-triangle HAF (viz., the J1-J2 
model) [^127) . This is a variant of the well-known J1-J2 
model on the infinite 2D square lattice, described below, 
in which one half of the NNN J2 bonds are removed. In 
the present paper we further the investigation of the J1-J2 
model by replacing the spin-i particles by particles with 
s = 1 and s = |. The 2D spin-i J1-J2 model has also 
been studied recently by other means [28l[29l[30ll31j . but 
we know of no other studies of the model for spins with 
spin quantum number s > i. 



\ 

Fig. 1: (colour online) J1-J2 model; - Ji; •^2; (a) Neel 

state, (b) spiral state, (c) striped state. 

2 The model 

The Hamiltonian of the J1-J2 model is written as 



if = Ji ^ s. 



[t,k] 



(1) 



where the operators = {s^,s^,sf) are the spin opera- 
tors on lattice site i, with sf = s{s + 1), and we consider 
here the cases 5 = 1 and s — ^. On the square lattice the 
sum over runs over all distinct NN bonds, but the 

sum over [i, k] runs only over one half of the distinct NNN 
bonds with equivalent bonds chosen in each square plaque- 
tte, as shown explicitly in Fig. [T] We shall be interested 
here only in the case of competing (or frustrating) antifer- 
romagnetic bonds Ji > and > 0, and henceforth for 
all of the results shown we set Ji = 1 . Clearly, the model 
may be described equivalently as a Heisenberg model on 
an anisotropic triangular lattice in which each triangular 
plaquette contains two NN Ji bonds and one NN bond. 
The model thus interpolates continuously between HAFs 
on a square lattice (J2 = 0) and on a triangular lattice 
( J2 = Ji). Similarly, when Ji = (or ^ 00 in our nor- 
malization with Ji = 1) the model reduces to uncoupled 
ID chains (along the chosen diagonals on the square lat- 
tice) . The case 3> 1 thus corresponds to weakly coupled 
ID chains, and hence the model also interpolates between 
ID and 2D scenarios. As well as the obvious theoretical 
richness of the model, there is also experimental interest 
since it is also believed to well describe such quasi-2D crys- 
talline materials as organic compounds containing BEDT- 
TTF [32 , for which with J2/J1 lies typically in the range 
from about 0.3 to about 1; and CS2CUCI4 [33], for which 
J2/J1 takes a value of about 6, thus making this material 
even quasi- ID. 

The J1-J2 model has only two gs phases in the classical 
case (corresponding to the limit where the spin quantum 
number s 00). For < 5 Ji the gs phase is Neel or- 
dered, as shown in Fig. [DJa), whereas for > ^Ji it 
has spiral order, as shown in Fig. (T^b), wherein the spin 
direction at lattice site points at an angle a^- = 



aa + (i + j)ach with ad 



cos 



i-Tk) ^ ^ - -^ci. The 



pitch angle (/)ci — cos ^{jjr) thus measures the deviation 

from Neel order, and it varies from zero for 2 Jj/ Ji < 1 to 
^TT as J2/J1 00, as shown later in Fig. [31 When = Ji 
we regain the classical 3-sublattice ordering on the trian- 
gular lattice with ad = |7r. The classical phase transition 
at J2 — ^ Ji is of continuous (second-order) type, with the 
gs energy and its derivative both continuous. 
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In the limit of large J'2/ Ji the above classical limit rep- 
resents a set of decoupled ID HAF chains (along the diag- 
onals of the square lattice) with a relative spin orientation 
between neighboring chains that approaches 90°. In fact, 
of course, there is complete degeneracy at the classical 
level in this limit between all states for which the relative 
ordering directions of spins on different HAF chains are 
arbitrary. Clearly the exact spin-^ limit should also be a 
set of decoupled HAF chains as given by the exact Bethe 
ansatz solution {34). However, one might expect that this 
degeneracy could be lifted by quantum fluctuations by the 
well-known phenomenon of order by disorder [3S]. Just 
such a phase is known to exist in the J1-J2 model [361137] 
for values of J2/ Ji ^ 0.6, where it is the so-called collinear 
striped phase in which, on the square lattice, spins along 
(say) the rows in Fig.[l]order ferromagnetically while spins 
along the columns and diagonals order antiferromagneti- 
cally, as shown in Fig. [IJc). We investigate the possibility 
below whether a stripe-ordered phase may be stabilized 
by quantum fluctuations at larger values of k for either of 
the cases s = 1 or s = |, in order to compare with the 
earlier s = ^ case for which we found that such a gs 
phase might exist for high enough values of the frustration 
parameter k, as discussed below. 

Thus, for the 5 = 5 case our own CCM calculations [26] 
provided strong evidence that the spiral phase becomes 
unstable at large values of the frustration parameter k. In 
view of that observation we also used the CCM for the s — 
^ case with the collinear stripe-ordered state as a model 
state. We found tentative evidence, based on the relative 
energies of the two states, for a second zero-temperature 
phase transition between the spiral and stripe-ordered states 
at a larger critical value of Kc^ ~ 1.8 ± 0.4, as well as firm 
evidence for a first phase transition between the Neel an- 
tiferromagnetic phase and the helical phase at a critical 
coupling = 0.80 ±0.01. 

The transition at k = for the s = case was found 
to be an interesting one. As in the classical (s — 00) case, 
the energy and its first derivative appeared to be contin- 
uous (within the errors inherent in our approximations), 
thus providing a typical scenario of a second-order phase 
transition, although a weakly first-order one could not be 
excluded since the gs energy did in fact show some defi- 
nite signs of a (weak) discontinuity in slope. Furthermore, 
the average on-site magnetization was seen to approach a 
value = 0.025 ± 0.025 very close to zero on both sides 
of the transition, but with a very sharp drop and hence a 
possible discontinuity in M on the spiral side of the tran- 
sition, as is often more typical of a first-order transition. 

A particular interest here is to compare and contrast 
the corresponding transition(s) between the s = and the 
s > ^ models. By contrast with the s = ^ we shall find 
below that for the cases with s = 1 and s = f the aver- 
age on-site magnetization at the analogous phase transi- 
tion between Neel-ordcrcd and spirally-ordered states ap- 
proaches smoothly the same nonzero value on both sides 
of the transition. Such continuous phase transitions where 
the order parameter does not vanish are well known in 
quantum magnetism. A prototypical example is the s — ^ 



anisotropic XXZ model with a Hamiltonian given by 

and which thus contains only NN anisotropic, antiferro- 
magnetic (J > 0) Heisenberg bonds. Classically (corre- 
sponding to the s ^ 00 limit) the model has a continuous 
phase transition a.t A = Ac = 1 between two different Neel 
antiferromagnetic phases, one aligned along the z-axis for 
A > 1, and the other along some arbitrary direction in 
the perpendicular xy-plane for — 1 < Z\ < 1. The s = ^ 
model may also be solved exactly on a ID chain by the 
Bethe ansatz technique [38]. It is found in this ID case 
that as the critical point is approached, A —i' A^ = 1, 
from either side, the average on-site (or staggered) mag- 
netization M Mc = 0. The approach to zero is of a 
quite nontrivial kind, via a function with an essential sin- 
gularity at Z\ = Ac- By contrast, for the same s = ^ XXZ 
model of Eq. (gj on a 2D square lattice, M ^ Mc ~ 0.31 
as Z\ —?> Z\c = 1. Thus, the type of phase transition we 
observe below for the spin-1 and s = f 2D interpolating 
square-triangle Heisenberg antiferromagnets is quite anal- 
ogous to the one a,t A = Ac = 1 in the s = ^ XXZ model 
on the 2D square lattice, but not to that of the same XXZ 
model on the ID chain. 

We now first briefly describe the main elements of the 
CCM below in Sec. O where we also discuss the approx- 
imation schemes used in practice for the s = ^ case and 
the s > I cases. Then in Sec. |4] we present our CCM re- 
sults based on using the Neel, spiral and striped states 
discussed above as model states (or starting states). We 
conclude in Sec. [5] with a discussion of the results. 

3 The coupled cluster method 

The CCM (see, e.g., Refs. [23,24,25 and references cited 
therein) is regarded as one of the most powerful and most 
versatile modern techniques in quantum many-body the- 
ory. It has been successfully applied to many quantum 
magnets (see Refs. l4tlMtlMt[37il391linillIlllgllMllili ) and 
references cited therein). The CCM is suitable for study- 
ing frustrated systems, for which the main alternative 
methods are often only of limited usefulness. For example, 
quantum Monte Carlo techniques are usually restricted by 
the sign problem for such systems, and the exact diagonal- 
ization method is limited in practice, especially for s > i, 
to such small lattices that it is often insensitive to the 
details of any subtle phase order present. 

The CCM method to solve the gs Schrodinger ket and 
bra equations, H\W) = E\^) and {W\H = respec- 
tively is now briefly outlined (and see Refs. 9 , 23 , 24, 25. BU1 
HO] for further details). The implementation of the CCM 
is initiated by the selection of a model state on top 
of which to incorporate later in a systematic fashion the 
multispin correlations contained in the exact ground states 
and The CCM employs the exponential ansatz, 
= e^\(l>) and {<P\ = (^|^e~'^. The creation correla- 
tion operator S is written as 5* = X]/=^o ^i^t with its 
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destruction counterpart as S* = 1 + "^j^qSjCj . The op- 
erators = {CyY, with Cq = 1, have the property 
that (<?|C/ = = Co\^) ; V/ 7^ 0. They form a com- 
plete set of multispin creation operators with respect to 
the model state |#). The calculation of the ket- and bra- 
state correlation coefficients (<S/,5/) is performed by re- 
quiring the gs energy expectation value H = {'F\H\'P) to 
be a minimum with respect to each of them. This results 
in a coupled set of equations {<P\CY e'^ He^ \<P) = and 
(<Z>|S'(e-^i/e^ - E)Cl\^) = ; V/ =^ 0, which we nor- 
mally solve by using parallel computing routines |45| for 
the correlation coefficients (5/,iS/) within specific trunca- 
tion schemes as outlined below. 

In order to treat each lattice site on an equal footing a 
mathematical rotation of the local spin axes on each lattice 
site is performed such that every spin of the model state 
aligns along its negative z-axis. As a result, our descrip- 
tion of the spins is given wholly in terms of these locally 
defined spin coordinate frames. The multispin creation op- 
erators may be expressed as = • • • , in terms 
of the locally defined spin-raising operators = sf + sf 
on lattice sites i. Upon solving for the multispin clus- 
ter correlation coefficients {Si,Si) as outlined above, the 
gs energy E may then be calculated from the relation 
E = {<1>\q~^ He^\<l>) ^ and the gs staggered magnetization 
M from the relation M = X]i=i in terms of 
the rotated spin coordinates. 

If a complete set of multispin configurations {/} with 
respect to the model state |^) is included in the calculation 
of the correlation operators S and 5, then the CCM for- 
malism becomes exact. However, it is necessary in practi- 
cal applications to use systematic approximation schemes 
to truncate them to some finite subset. For the s = 5 case, 
the localised LSUBn scheme is commonly employed, as in 
our earlier paper on the 5 = 5 version of the present model 
[26] . as well as in our other previous work [9 > 25, 39, 40 > 43 . 
Under this truncation scheme all possible multi-spin-fiip 
correlations over different locales on the lattice defined by 
n or fewer contiguous lattice sites are retained. A cluster 
is defined as having n contiguous sites if every one of the 
n sites is adjacent (as a nearest neighbour) to at least one 
other. Clearly this definition, however, depends on how we 
choose the geometry of the lattice among various topolog- 
ically equivalent possibilities that may exist. For example, 
the current model may be construed as referring to sites 
on a square lattice, as shown in Fig. [TJ In this case the 
bonds, for example, join NNN sites (which, by definition, 
are thus not adjacent). Alternatively, the model may be 
equivalently construed as referring to sites on a triangular 
lattice, in which case both Ji and J'2 bonds join NN (and 
hence adjacent) sites. In all of the results presented here 
we consider the model to be defined on a triangular lattice 
in making CCM approximations. 

However, we note that the number of fundamental 
LSUBn configurations for s > ^ becomes appreciably 
higher than for s = i, since each spin on each site i can 
now be raised up to 2s times by the spin-raising oper- 
ator . Thus, for the s > ^ models it is more practi- 



Table 1: Number of fundamental CCM configurations 
{Nf) for the SUBn-n {n = {2, 3, 4, 5, 6, 7}) scheme for the 
Neel and striped model states, for the J1-J2 model defined 
on a triangular lattice, for the spin-1 and spin-| cases. 



Method 


s = 


1 


s = 


3 
2 


Nf 




Nf 




SUBn-n 


striped 


spiral 


striped 


spiral 


SUB2-2 


2 


4 


2 


4 


SUB3-3 


4 


26 


4 


27 


SUB4-4 


60 


189 


60 


211 


SUB5-5 


175 


1578 


175 


1908 


SUB6-6 


2996 


14084 


3622 


18501 


SUB7-7 


11778 


131473 


13320 


188326 



cal, but equally systematic, to use the alternative SUBn- 
m scheme, in which all correlations involving up to n 
spin flips spanning a range of no more than m contigu- 
ous lattice sites are retained [TSlfTBlfTB II 2 5 [HT] . We then set 
TO = n, and hence employ the so-called SUBn-n scheme. 
More generally, the LSUBto scheme is thus equivalent to 
the SUBti-to scheme for n — 2sm for particles of spin 
s. For s — ^, LSUBn = SUBn-n; whereas for s > i, 
LSUBn = SUB2sn-n. The numbers of such fundamen- 
tal configurations (viz., those that are distinct under the 
symmetries of the Hamiltonian and of the model state |#)) 
that are retained for the Neel and striped model states of 
the current s = 1 and = | models at various SUBn-n 
levels, defined with respect to an underlying triangular- 
lattice geometry, are shown in Table [T] 

Although we never need to perform any finite-size scal- 
ing, since all CCM approximations are automatically per- 
formed from the outset in the N 00 limit, where N is 
the total number of lattice sites, we do need as a last step 
to extrapolate to the n ^ 00 limit in the truncation index 
n. We use here the well-tested [151116112611401141] empirical 
scaling laws 

E/N = ao + aiij-'^ + a2n-^ , (3) 

M = bo + bin-^ + b2n-^ , (4) 

exactly as we did previously for the corresponding s = ^ 
model pS], for the gs energy per spin E/N and the gs 
staggered magnetization M, respectively. 

4 Results 

The results of the CCM calculations are reported here 
for the spin-1 and spin-| J1-J2 model Hamiltonian of Eq. 
(HI, using the Neel, spiral and striped states shown in Fig, 
l(a)-(c) as CCM model states, and with the SUBn-n ap- 
proximation scheme defined with respect to an underlying 
triangular-lattice geometry. We set the parameter Ji — 1. 
Our available computational power at present is such that 
we can perform SUBn-n calculations for the spiral model 
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state (viz., the state that requires the highest number of 
fundamental configurations (Nf) for a given SUB71-71 trun- 
cation index n) only for values n < 7 for both the s — 1 
and s = I cases. We thus present results for each of the 
Neel, striped and spiral states only up to the SUB7-7 level, 
for the sake of consistency in our extrapolations to the 
n — 00 limit . 

We note that, as has been well documented in the 
past [3^, the LSUBn (or SUBn-n) data for both the gs 
energy per spin E/N and the average on-site magneti- 
zation M converge differently for even-n sequences and 
odd-n sequences, similar to what is frequently observed 
in perturbation theory [4 7) . Since, as a general rule, it is 
desirable to have at least (n + l) data points to fit to any 
fitting formula that contains n unknown parameters, we 
prefer to have at least 4 results to fit to Eqs. ^ and dH). 
Both the available odd and even series of our SUBrt-n data 
violate this desirable rule. However, our results (for both 
sets n = {2, 4, 6} and n = {3, 5, 7}) for the s = ^ case are 
consistent with those using the larger LSUBto sequences 
available in this case. This gives us confidence in both the 
accuracy of our results and the robustness of our extrapo- 
lation schemes. Hence, for most of our extrapolated results 
below we use the even SUBn-n sequence with n — {2, 4, 6} 
and the odd SUBn-n sequence with n = {3, 5, 7}. 

Firstly, the results obtained using the spiral model 
state are reported. For this state we first perform CCM 
calculations with the pitch angle ^ as a free parameter. At 
each separate level of approximation we then choose the 
angle 4> = 4>SVBn-n that minimizes the energy EsuBn-n{4>)- 

Classically we have a second-order phase transition 
from Neel order (for k < Kci) to helical order (for k > Kd), 
where k = J2I Ji, at a value Kd = 0.5. By contrast, our 
CCM results presented below show that there is a shift 
of this critical point to a value Kc ~ 0.615±0.010 in the 
spin-1 quantum case and Kc ~ 0.575 ±0.005 for the spin-| 
quantum case, first indications of which are seen in Figs. 

and 121 In both cases this is a second-order phase transi- 
tion from Neel-ordered to helically-ordered states. Thus, 
for example, curves such as those shown in Fig. [5] show 
that the Neel state (0 = 0) gives the minimum gs energy 
for all values of k < Kc, where Kc depends on the level of 
SUBn-n approximation used, as we also observe in Fig. 
[3l By contrast, for values of k > Kc the minimum in the 
energy is found to occur at a value (/i 7^ 0. If we consider 
the pitch angle (j) itself as an order parameter (i.e., = 
for Neel order and ^ for spiral order) a typical sce- 
nario for a phase transition would be the appearance of 
a two-minimum structure for the gs energy for values of 
K > Kc, exactly as observed in Fig. [5] for both the spin- 
1 and spin- 1 models in the SUB4-4 approximation. Very 
similar curves occur for other SUBn-n approximations. 

We note that the crossover from one minimum (0 = 0, 
Neel) solution to the other (0 ^ 0, spiral) appears to be 
quite smooth at this point (and see Figs. [2] and [3]). Thus, 
for example, the spiral pitch angle 4> appears to change 
quite continuously from a value of zero for k < on the 
Neel side of the transition to a nonzero value for k > Kc 
on the spiral-phase side. For example, at the SUB6-6 level 



we find Kc ~ 0.613 for the spin-1 case, and Kc ~ 0.574 
for the spin-spin-| case. We also note from Fig. [3] that 
as J2 — > 00 the spiral angle approaches the limiting 
value ^TT considerably slower for the spin-1 and spin-| 
cases than it does the spin-i case we investigated earlier 
(and see Fig. 3 in Ref. [Ml)- This is a first indication that 
there is less freedom for the existence of a stable collinear 
(striped) state at higher values of k for the higher spin 
(s > ^) models than for the 5 = 5 model. We return to 
this point later. 

Figure [2] shows the ground-state energy per spin versus 
the spiral angle 0, using the SUB4-4 approximation of the 
CCM with the spiral model state, for some illustrative 
values of J2. Similarly Fig.[3]shows the angle (psvBn-n that 
minimizes the energy £'suBn-n(0)- Our previous study of 
the quantum spin-^ case in the same model |26| found 
that there is a first quantum critical point at w 0.80 
at which a weakly first-order, or possibly second-order, 
phase transition occurs between states that exhibit Neel 
order and helical order. We see now that increasing the 
spin quantum number s thus brings the quantum critical 
point Kc closer to the classical critical point k^ — 0.5 for 
the phase transition from Neel order to helical order, as 
expected. 

We observe from Fig. [2] that for certain values of 
(or, equivalently, k) CCM solutions at a given SUBn-n 
level of approximation (viz., SUB4-4 in Fig. exist only 
for certain ranges of the spiral angle (j). For example, for 
the pure square-lattice HAF {k = 0) the CCM SUB4- 
4 solution based on a spiral model state only exists for 
< ^ 0.177r for the spin-1 model and < < O.Wtt 
for the spin- 1 model. In this case, where the Neel solu- 
tion is the stable ground state, if we attempt to move too 
far away from Neel collinearity the CCM equations them- 
selves become "unstable" and simply do not have a real 
solution. Similarly, we see from Fig.[2]that for k = 1.5 the 
CCM SUB4-4 solution exists only for 0.257r < (p < O.Btt 
for the spin-1 model and for 0.277r < (f> < 0.5tt for the 
spin- 1 model. In this case the stable ground state is a spi- 
ral phase, and now if we attempt to move too close to Neel 
collinearity the real solution terminates. 

Such terminations of CCM solutions are common |25) . 
A termination point usually arises because the solutions 
to the CCM equations become complex at this point, be- 
yond which there exist two branches of entirely unphysical 
complex conjugate solutions |25| . In the region where the 
solution reflecting the true physical solution is real there 
actually also exists another (unstable) real solution. How- 
ever, only the (shown) upper branch of these two solutions 
reflects the true (stable) physical ground state, whereas 
the lower branch does not. The physical branch is usually 
easily identifled in practice as the one which becomes ex- 
act in some known (e.g., perturbative) limit. This physical 
branch then meets (with infinite slope, as seen in Fig. ^ 
the corresponding unphysical branch at some termination 
point beyond which no real solutions exist. The SUBn-n 
termination points are themselves also reflections of the 
quantum phase transitions in the real system, and may 
be used to estimate the position of the phase boundary 
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pS] . although we do not do so here since we have more 
accurate criteria discussed below. 

Figures|l]and[S]show the CCM results for the gs energy 
and average gs on-site magnetization, respectively, where 
the spiral state has been used as the model state. The 
gs energy (in Fig. 2]) shows no sign of a discontinuity in 
slope at the critical values Kc discussed above, and this is 
an indication of a second-order transition from the Neel 
phase to the helical phase. This is in contrast with the 
spin-i case, where the gs energy shows definite signs of a 
(weak) discontinuity in slope at the first critical value 

m- 

The gs magnetic order parameter M in Fig. [S] shows 
much clearer evidence of a phase transition at the corre- 
sponding Kc values previously observed in Fig.|31 Thus, we 
see that for the spin-1 case the sharp minimum in the ex- 
trapolated magnetic order parameter occurs at Kc ~ 0.613 
(with Mc = 0.6367) using n = {2,4,6} for the extrap- 
olation of M, and at ~ 0.606 (with Mc = 0.5978) 
using n = {3, 5, 7}; whereas for the spin-| case, the corre- 
sponding values are Kc ~ 0.574 (with Mc = 1.1134) using 
n = {2, 4, 6} for the extrapolation of M and at Kc « 0.571 
(with Mc — 1.0766) using n — {3,5,7}. We also present 
other independent estimates for Kc below. 

By contrast, for the spin-i case j^S] the extrapolated 
value of M showed clearly its steep drop toward a value 
very close to zero at a corresponding value Kc ~ 0.80, 
which gave the best CCM estimate of the phase-transition 
point for that case. In the spin-i case the magnetization 
seemed to approach continuously a value M = 0.025 ± 
0.025 from the Neel side {k < Kc) whereas from the spiral 
side (k > Kc) there appeared to be a discontinuous jump 
in M as K — >■ Kc. The transition at k < Kc thus appeared 
to be (very) weakly first order but it was not possible to 
exclude it being second order since the possibility of a 
continuous but very steep drop to zero of the on-site mag- 
netization as K ^ Kc from the spiral side of the transition 
could not be entirely ruled out. No evidence at all was 
found for any intermediate phase between the quasiclassi- 
cal Neel and spiral phases, just as for the higher-spin cases 
considered here. However, Fig. [5] here shows no evidence 
at all for a finite jump in M as k — > Kc from either side 
of the transition, and hence the evidence from the order 
parameter is that the transition from Neel order to spiral 
order for both the spin-1 and spin-| cases is of second- 
order type. 

Table [5] shows the critical values KSUBn-n which 
the transition between the Neel and spiral phases occurs 
in the various SUBn-n approximations shown in Fig. |3l 
In the past we have found that a simple linear extrap- 
olation scheme HHSlIll], KSUBn-™ = ao + ain"!, yields 
a good fit to such critical points. This seems to be the 
case here too, just as for the spin-i case [35]. The fact 
that the two corresponding "SUBoo" estimates from the 
SUBn-n data in Table [5] based on the even-n and odd- 
n SUBn-n sequences differ slightly from one another is a 
reflection of the errors inherent in our extrapolation pro- 
cedures. Similar estimates based on an alternative extrap- 
olation scheme, kLSUBji _ _|_ l)^ri~'^, are also shown in 



Table 2: The critical value Kc = at which the 

transition between the Neel phase (0 = 0) and the spiral 
phase {(j) 7^ 0) occurs in various SUBn-n approximations, 
using the CCM with the (Neel or) spiral state as model 
state, for the J1-J2 model. Results are shown for both the 
spin-1 and spin-| cases. 



Method 


s = 1 


3 






Hie 


SUB2-2 


0.597 


0.563 


SUB4-4 


0.610 


0.571 


SUB6-6 


0.613 


0.574 


SUBoo " 


0.617 


0.581 


SUBoo 


0.616 


0.577 


SUB3-3 


0.577 


0.554 


SUB5-5 


0.597 


0.566 


SUB7-7 


0.607 


0.571 


SUBoo " 


0.636 


0.583 


SUBoo 


0.619 


0.577 


" Based on 


1/n : n — 


{2,4,6} 


^ Based on 


: n = 


= {2,4,6} 


Based on 


1/n : n — 


{3,5,7} 


Based on 


1/n^ : n = 


= {3,5,7} 



Table [D The difference between all of these estimates is 
thus also a rough indication of our real error bars on Ke- 
lt is gratifying to note that all of the estimates for 
Kc from the extrapolations of our computed results for 
^SUBn-n excellent agreement with those obtained 

from the extrapolated results for the order parameter M 
discussed above. By putting all of these results together, 
our final estimates for the critical point for the transition 
between the Ncel-ordered and the spirally-ordered phases 
are Kc = 0.615 ± 0.010 for the spin-1 model and Kc = 
0.575 ± 0.005 for the spin-§ model. 

We conclude our discussion of the Neel and spiral phases 
by presenting detailed results for the two spin cases for the 
two special limits of the model, namely the pure isotropic 
HAF on the square and triangular lattices. Thus, Table |3] 
shows the results for the ground-state energy per spin and 
magnetic order parameter (i.e., the average on-site mag- 
netization) for the spin-1 and spin-| J1-J2 HAF model on 
the square lattice (J2 = or k — 0) and on the trian- 
gular lattice (J2 = Ji or K = 1), using the spiral model 
state. Our CCM results are presented in various SUBn-n 
approximations (with 2 < n < 7) based on the triangular 
lattice geometry using the spiral model state, with (j) = Q 
for the square lattice and </) = ^ for the triangular lattice. 
The extrapolated results (n 00) using Eqs. ^ and ^ 
with n = {2, 4, 6} and n = {3, 5, 7} are also presented. For 
comparison we also show the results obtained for the spin- 
1 model on the square lattice (i.e., k = 0) using spin- wave 
theory (SWT) [48j, a linked-cluster series expansion (SE) 
method ^49j, and previous CCM SUBn-n (n 00) results 
based on the model construed as referring to sites on a 
square lattice T3]. Our present results are seen both to 
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Table 3: Ground-state energy per spin and magnetic order parameter (i.e., the average on-site magnetization) for the 
spin-1 and spin-| HAFs on the square and triangular lattices. We show COM results obtained for the J1-J2 model 
with Ji > 0, using the spiral model state in various SUBn-n approximations defined on the triangular lattice geometry, 
for the two cases k = J'lj J\ = (square lattice HAF, </) = 0) and k = 1 (triangular lattice HAF, = -1). 



iVietnoa 


E/N 


M 


E/N 


M 


E/N 


M 


E/N 


M 




square 


(^ = 0) 


triangular (k = 1) 


square (k 


= 0) 


triangular 




SUB2-2 


-2.29504 


0.9100 


-1.77400 


0.9069 


-4.94393 


1.4043 


-3.80006 


1.3938 


SUB3-3 


-2.29763 


0.9059 


-1.80101 


0.8791 


-4.94836 


1.3990 


-3.83393 


1.3672 


SUB4-4 


-2.31998 


0.8702 


-1.82231 


0.8405 


-4.97694 


1.3638 


-3.86025 


1.3287 


SUB5-5 


-2.32049 


0.8682 


-1.82623 


0.8294 


-4.97789 


1.3611 


-3.86498 


1.3170 


SUB6-6 


-2.32507 


0.8510 


-1.83135 


0.8096 


-4.98305 


1.3452 


-3.87059 


1.2980 


SUB7-7 


-2.32535 


0.8492 


-1.83288 


0.8006 


-4.98344 


1.3430 


-3.87191 


1.2904 



Extrapolations 



SUBoo " 


-2.32924 


0.8038 


-1.83860 


0.7345 


-4.98793 


1.3001 


-3.87869 


1.2233 


SUBoo '' 


-2.32975 


0.7938 


-1.83968 


0.7086 


-4.98803 


1.2933 


-3.87839 


1.2107 



CCM " 
SWT 
SE = 



-2.3291 
-2.3282 
-2.3279(2) 



0.8067 
0.8043 
0.8039(4) 



" Based on n = {2, 4, 6} 
Based on n = {3, 5, 7} 

CCM (SUBoo for square lattice, based on n = {2, 4, 6}) in the natural square-lattice geometry [13] 
SWT (Spin-wave theory) for square lattice [48] 
" SE (Series Expansion) for square lattice [15] 



be robust and internally consistent, by comparison of the 
independent extrapolations of the SUBn-n data using the 
even-n and odd-n data sets, and to agree very well with 
the best alternative results available for the spin-1 model 
on the square lattice. Such comparisons give us confidence 
that our results are likely to be similarly accurate over the 
entire range of values of the frustration parameter k. 

We turn finally to our CCM results based on the collinear 
striped AFM state as the choice for the CCM gs model 
state \<1>). The SUBn-n configurations are again defined 
with respect to the triangular lattice geometry, exactly 
as before. The numbers of fundamental configuration Nf 
in each of the SUBn-n approximations used are given in 
Table [1] Results for the gs energy and magnetic order pa- 
rameter based on the striped phase are shown in Figs. [6] 
and [7] respectively. We see from Fig. [S] that some of the 
SUBn-n solutions based on the striped state for both the 
s = 1 and s = § cases show a clear termination point Kt 
of the sort discussed previously, such that for k < Kt no 
real solution for the striped phase exists. 

For the spin-1 model the large-K limit of the extrapo- 
lated SUBn-n energy per spin results of E/N — —1.3897 
from Fig.[ni;a) using n = {2,4,6} and E/N = -1.3936 
from Fig. ^h) using n — {3, 5, 7} agree well with the 
known ID chain result of E/N ~ —1.4015 obtained from 
a density-matrix renormalization group analysis |50| and 
our previous CCM result 9 , just as in Fig. [4]Ja) and (b) 
for the spiral phase. Similarly, for the spin-| case, the 
large-K limit of the extrapolated SUBn-n results for the 



energy per spin of E/N = — 2.8205J2 from Fig.[ni^c) using 
n = {2,4,6} and E/N = -2.8243 J;^ from Fig. [Bjjd) using 
n = {3, 5, 7}, with almost identical results again obtained 
from Fig. \^c) and (d). Unlike for their spin-i counter- 
part, however, the striped phase is never a stable gs state 
for either the spin-1 or spin-| models, because their ener- 
gies always lie higher than those of the spiral state for all 
values of J2, as shown in Fig. [5] Hence for the s = 1 and 
s = I cases, there is only one quantum critical point Kc, 
at which the Neel phase is driven to the helical phase. 



5 Discussion and conclusions 

In an earlier paper |2S] we used the CCM to study the 
effect of quantum fluctuations on the zero-temperature gs 
phase diagram of a frustrated spin-i interpolating square- 
triangle antiferromagnetic model. This is the so-called Ji- 
J2 model, defined on an anisotropic 2D lattice, as shown 
in Fig. [TJ In the current paper we have extended the anal- 
ysis to consider spin-1 and spin-| versions of the same 
model. As before we have studied the case where the NN 
Ji bonds are antiferromagnetic ( Ji > 0) and the compet- 
ing J2 = kJi bonds have a strength k that varies from 
K = (corresponding to the HAF on the square lattice) 
to K — ?► 00 (corresponding to a set of decoupled ID HAF 
chains) , with the HAF on the triangular lattice as another 
special 1, in between the two extremes. The re- 

sults of the K ~ limit of the present model (and see 
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Table [2]) for the s = 1 case are comparable with those ob- 
tained from the SWT and SE techniques [48ll49] which are 
among the best alternative numerical method to the CCM 
for highly frustrated spin-lattice models like the present 
J1-J2 model. 

For the spin-1 model we find that the phase transition 
between the Neel antiferromagnetic phase and the spiral 
phase occurs at the value Kc = 0.615 ± 0.010, whereas 
for the spin- 1 model we find that the phase transition 
occurs at Kc — 0.575 ± 0.005. From the continuous and 
smooth behaviour of the energies of the two phases it ap- 
pears that the transition is second-order, as in the classical 
case. However, on neither side of the transition at Kc does 
the order parameter M (i.e., the average on-site magneti- 
zation) go to zero for either of the two higher spins con- 
sidered here. On the other hand, unlike in the spin-^ case, 
in neither of the higher-spin models does there appear to 
be any discontinuity in M at the transition. All of the 
indications are thus that the transition between the Neel 
antiferromagnetic and the spiral phases is of continuous 
(second-order) type for both cases s = 1 and s = |, in 
contrast to the spin-^ case where the order parameter M 
appeared to show a discontinuous jump at the transition, 
which was found to be a weakly first-oder one (although 
it could not be entirely excluded on the available evidence 
that the transition might be a second-order one). 

We have observed that as the quantum spin number 
s is increased, the position of the quantum critical point 
at Kc between the phases with Neel and spiral order is 
brought closer to the classical (s — 00) value, KcI = 0.5, 
as expected. In contrast with the s — case where there 
is a second quantum critical point for the phase transition 
from the helical phase to a collinear stripe-ordered phase, 
we find no evidence at all for such a further transition for 
either of the cases s = 1 or s = | . 

We note that the spin-1 HAF on the (undistorted) tri- 
angular lattice (viz., our limiting case k = 1) has itself 
been the subject of much recent interest from both the 
theoretical and experimental viewpoints. From the exper- 
imental side spin-1 models on the triangular lattice are 
believed to underlie the properties of such materials as 
NiGa2S4 [SU and Ba3NiSb209 [52 . In both materials the 
Ni2+ ions form in weakly coupled 2D triangular lattice 
layers. Thus, for example, thermodynamic and neutron 
scattering measurements on NiGa2S4 show conclusive ev- 
idence that the inherent geometric frustration of the trian- 
gular lattice stabilises a low-temperature spin-disordered 
state, which was proposed as being consistent with a spin- 
liquid phase 51]. Other candidates for spin-1 quantum 
spin-liquid phases have more recently been proposed from 
an experimental study of the high-pressure sequence of 
structural phases in the material Ba3NiSb209 [55]. 

Whereas quantum fluctuations are certainly intrinsi- 
cally greatest for spin-lattice systems with the lowest spin 
value s = i, as we have also found here, such effects can 
also be enhanced for the s > 5 cases by the addition to 
the pure (bilinear) Heisenberg interaction with NN terms 
only of terms such as a NN biquadratic interaction or 
other higher-order exchange terms. It is precisely by the 



addition of terms like this that unusual quantum ground 
states, such as ones with quadrupolar (or spin-nematic) 
order have been predicted theoretically to be stabilised 
for the spin-1 HAF on the triangular lattice [551I5HI551 
[55] . It is argued that such a state can account for the 
observed low-temperature thermodynamics in the spin-1, 
quasi-2D, antiferromagnetic material NiGa2S4, although 
at the lowest temperatures the observed order is that of 
an (incommensurate) spiral phase. In a very recent pa- 
per (that appeared only after submission of this paper) 
|57) it is also argued that the quantum spin-liquid phases 
presumed to have been seen in recent experiments [55] 
in the layered material Ba3NiSb209, may be explained 
microscopically as emanating from a spin-1 HAF on the 
triangular lattice with both NN and NNN isotropic anti- 
ferromagnetic Heisenberg couplings. 

Such other (e.g., spin-1) models on the triangular lat- 
tice as those described above, involving either isotropic 
bilinear and biquadratic couplings or both NN and NNN 
bilinear Heisenberg couplings, could also be investigated 
via the CCM, and it would surely be interesting to do so. 
While such additional terms in the Hamiltonian present 
no additional obstacles to the use of the method at all, 
the choice of which model states to use always enters at 
the outset. It is certainly true that most calculations on 
spin systems employing the CCM, including those in the 
present paper, employ model states built by independent- 
spin product states for which the choice of state for the 
spin on each site is formally independent of the choice 
of all others. Often for these independent-spin product 
model states the use of collinear states, such as the Neel 
or striped states considered here, is possible, where all 
spins are aligned parallel or antiparallel to one axis. How- 
ever, as we have seen, noncoUinear (e.g. spiral) model 
states can sometimes be favourable for certain values of 
the frustration. In either case multispin correlations are 
then included systematically on top of the independent- 
spin product model states. As we have seen here, the CCM 
for such independent-spin product model states may then 
be applied to high orders by using a computational imple- 
mentation used here and described more fully elsewhere 
(see, e.g., Refs. [391191145] and references cited therein). In 
particular, it may be applied to lattices of complex crys- 
tallographic symmetry. Furthermore, as seen here, it is not 
constrained to systems with spin quantum number s = ^. 

When the system under consideration may have more 
exotic ground states with less conventional ordering than 
the (often essentially quasiclassical) independent-spin prod- 
uct states described above, the CCM may still be very 
profitably employed. Even the use of such independent- 
spin product states can still give very precise phase bound- 
aries for when such states give way to more exotic states. 
A good example among many to date is the well-studied 
frustrated spin-J J1-J2 model on the square lattice dis- 
cussed in Sec. [J for which the phase boundaries of the 
non-classical paramagnetic state (that has no magnetic 
LRO) have been estimated very accurately using the CCM 
in Refs. [551157] . A more recent example is provided by a 
CCM calculation [5F| of the frustrated spin-^ J1-J2-J3 
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model on the honeycomb lattice which incorporates NN 
bonds (Ji) and NNN bonds (J2) as in the J1-J2 model, 
but now also includes next-next-nearest-neighbour bonds 
(J3). For the case J3 = J2 an intermediate paramagnetic 
phase was accurately located between collinear antiferroi- 
magnetic states of quasiclassical Neel and striped order. 
By calculating with such model states the plaquette sus- 
ceptibility, the authors gave precise values not only of the 
phase boundaries of this intermediate state, but also gave 
clear evidence that it had plaquette valence-bond crys- 
talline ordering. 

It is also worth noting that the CCM can deal directly 
with more complex model states, such as those involv- 
ing valence-bond crystal (VBC) order. Thus, for exam- 
ple, non-classical VBC ordering has been considered us- 
ing the CCM by employing directly valence-bond model 
states, i.e. two- or multi-spin singlet product states [59 . 
A drawback of this approach is that it involves the direct 
use of products of localized states (e.g., two-spin dimers 
or multi-spin plaquettes) in the model state. Hence, this 
approach requires that a new matrix-operator formalism 
be created for each new problem. Also, the Hamiltonian 
and CCM ket- and bra-state operators must be written 
in terms of this new matrix algebra. The CCM equations 
may be derived and solved once the commutation rela- 
tionships between the operators have been established. 
Although formally straightforward, this process can be 
tedious and time-consuming. Furthermore, the existing 
high-order CCM formalism and codes also need to be 
amended extensively for each separate model considered. 

More recently a quite different CCM approach has 
been advocated for dealing with such VBC states [60 . 
It starts directly from collinear independent- spin product 
model states, and shows how one may form exact local 
dimer or plaquette ground states within the CCM frame- 
work. This approach has the huge advantages of being 
conceptually simple and thus also of being easy to imple- 
ment. Furthermore, one may then use directly the existing 
high-order CCM formalism, computer codes, and extrap- 
olation schemes used and described here and in that ref- 
erences cited. To date the method has been applied with 
excellent results to the spin-i J1-J2 model for the lin- 
ear chain, the spin-i Shastry-Sutherland model on the 2D 
square lattice f&V, and the so-called spin-^ J- J' HAF on 
the 2D CAVO lattice that is appropriate to the magnetic 
material CaV409. It is a one- fifth depleted square lattice, 
and the model on this lattice comprises two nonequivalent 
antiferromagnetic NN bonds of strength J and J'. The J 
bonds connect sites on the NN four-spin square plaquettes 
while the J' (dimer) bonds connect NN sites belonging to 
neighbouring square plaquettes. 

In conclusion, it will be of interest to use the CCM for 
the other spin-1 models discussed above on the triangular 
lattice that are believed to be relevant to such quasi-2D 
materials as NiGa2S4 and Ba3NiSb209, for both of which 
considerable experimental data exist. We hope to be able 
to perform and report ourselves on such calculations at a 
later date. 
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Fig. 2: (colour online) Ground-state energy per spin of the spin-1 and spin-| J1-J2 Hamiltonian of Eq. ([T]) with Ji = 1, 
using the SUB4-4 approximation of the CCM with the spiral model state, versus the spiral angle </>. For the case of 
s = 1, for J2 ^ 0.610 the only minimum is at = (Neel order), whereas for J2 ^ 0.610 a secondary minimum occurs 
at = 0SUB4-4 7^ 0, which is also a global minimum, thus illustrating the typical scenario of a second-order phase 
transition. Similarly, for the case of s = |, for < 0.571 the only minimum is at (/> = (Neel order), whereas for 
J2 > 0.571 a secondary minimum occurs at (p = (/>suB4-4 7^ 0. 
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(a) s = 1 (b) s = 1 
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Fig. 3: (colour online) The angle ^suBn-n that minimizes the energy -EsuBn-n(0) of the spin-1 and spin-| J1-J2 
Hamiltonian of Eq. ([IJ with Ji = 1, in the SUBn-n approximations with n — {2,4,6} and n = {3,5,7}, using the 
spiral model state, versus J2. The corresponding classical result 0ci is shown for comparison. We find in the SUBn-n 
quantum case a second-order phase transition (e.g., at the SUB6-6 level, at w 0.613 for the s = 1 case and at 
J2 « 0.574 for the s = | case), where 0suBn-n changes continuously from zero below the transition point (Neel phase) 
to a nonzero value above it (helical phase). The classical case has a second-order phase transition at J2 = 0.5. 
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(c) s = I id) 3 = 1 

Fig. 4: (colour online) Ground-state energy per spin versus J2 for the Neel and spiral phases of the spin-1 and spin-| 
J1-J2 Hamiltonian of Eq. ([1]) with Ji = 1. The CCM results using the spiral model state are shown for various SUBn-rt 
approximations (n = {2, 4, 6}) and (n = {3, 5, 7}) with the spiral angle (j) = 0suBn-n that minimizes i?suBn-ri(0)- We 
also show the n — 00 extrapolated results from using Eq. ([3]). 
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(c) s = I id) 3 = 1 

Fig. 5: (colour online) Ground-state magnetic order parameter (i.e., the average on-site magnetization) versus J2 for 
the Neel and spiral phases of the spin-1 and spin-| J1-J2 Hamiltonian of Eq. ([1]) with Ji ~ 1. The CCM results using 
the spiral model state are shown for various SUBn-n approximations [n — {2, 4, 6} and {n = {3, 5, 6}) with the spiral 
angle (p = 4>svBn-n that minimizes i?suBn-n(0)- We also show the n ^ 00 extrapolated results from using Eq. (g]). 
The sharp minimum in the extrapolated magnetic order parameter is at J2 = 0.613 (M — 0.6367) using n = {2, 4, 6} 
and J2 = 0.606 (M = 0.5978) using n — {3, 5, 7} for the spin-1 case, whereas for the spin-| case, the corresponding 
values are = 0.574 (M = 1.1134) using n = {2,4,6} and = 0.571 (M = 1.0766) using n = {3,5,7}. 
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Fig. 6: (colour online) Ground-state energy per spin versus J'2 for the stripe-ordered phase of the spin-1 and spin-| 
J1-J2 Hamiltonian of Eq. ^ with Ji — 1. The CCM results using the striped model state are shown for various 
SUBn-n approximations {n = {2, 4, 6}) and (n = {3, 5, 7}). We also show the n — >■ 00 extrapolated results from using 
Eq. ©. 
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Fig. 7: (colour online) Ground-state magnetic order parameter (i.e., the average on-site magnetization) versus J 2 for 
the stripe-ordered phase of the spin-1 and spin-| Ji- J2 Hamiltonian of Eq. ([ij with Ji — 1. The CCM results using 
the striped model state are shown for various SUBrt-n approximations (n — {2, 4, 6} and n — {3, 5, 7}). We also show 
the n 00 extrapolated results from using Eq. 
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Fig. 8: (colour online) Difference between the ground-state energies per spin (e = E/N) of the spiral and striped 
phases {Ae = e'^P"'^' - e''*"?'''^) versus for the spin-1 and spin-| Ji-J^ Hamihonian of Eq. (H]) with Ji = 1. The CCM 
results for the energy difference using both the striped and spiral model states for various SUBn-n approximations 
{n = {2, 4, 6}) and (n = {3, 5, 7}) are shown. We also show the n — oo extrapolated results from using Eq. ^ for the 
two phases separately. 



